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1 .  INTRODUCTION 

The  problems  of  testing  the  hypotheses  on  the  structures 
of  the  covariance  matrices  of  the  real  multivariate  normal 
populations  received  considerable  attention  in  recent  years. 

But,  not  much  work  was  done  on  the  covariance  structures  of 
the  complex  multivariate  normal  populations.  Invest igat ions 
on  covariance  structures  of  the  complex  multivariate  normal 
populations  have  important  applications  in  the  area  of  inference 
on  multiple  time  series  since  certain  estimates  of  the  spectral 
density  matrix  of  the  multivariate  stationary  Gaussian  time 
series  are  approximately  distributed  as  the  complex  Wishart 
matrix.  The  object  of  this  paper  is  to  investigate  the  null 
and  nonnull  asymptotic  distributions  of  the  likelihood  ratio 
statistics  for  testing  the  hypotheses  on  the  covariance  struc¬ 
tures  of  the  complex  multivariate  normal  populations.  For 
some  discussions  on  the  applications  of  the  complex  multivariate 
distributions,  the  reader  is  referred  to  Brillinger  (1974), 
Hannan  (1970)  and  Krishnaiah  (1976). 

In  Section  3  of  this  paper,  we  derived  the  asymptotic  null 
distribution  of  the  likelihood  ratio  statistic  for  multiple 
independence  whereas  Section  4  is  devoted  to  the  corresponding 
distribution  in  the  nonnull  case  under  certain  alternatives. 

An  expression  is  derived  in  Section  5  for  the  asymptotic  null 
distribution  of  the  likelihood  ratio  statistic  for  multiple 
homogeneity  of  the  covariance  matrices  of  the  complex  multi¬ 
variate  normal  populations.  In  Section  6,  we  derived  an  ex- 
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pression  for  the  asymptotic  nonnull  distribution  of  the  like¬ 
lihood  ratio  test  statistic  for  homogeneity  of  the  covariance 
matrices  of  two  complex  multivariate  normal  populations  under 
certain  alternatives.  The  hypotheses  considered  in  Sections 
3-6  arise  in  studying  certain  linear  structures  of  the  co- 
variance  matrices.  For  a  discussion  of  these  problems  in  th* 
real  case,  the  reader  is  referred  to  Krishnaiah  and  Lee  (1974). 
The  expressions  obtained  in  this  paper  are  in  terms  of  the 
beta  series.  In  the  null  cases,  it  is  found  that  the  accuracy 
of  the  approximations  based  on  the  first  terms  of  the  asymptotic 
series  are  sufficient  for  many  practical  purposes.  Krishnaiah, 
Lee  and  Chang  (1976)  approximated  the  null  distributions  of 
certain  powers  of  the  likelihood  ratio  test  statistics  for 
multiple  independence  and  multiple  homogeneity  of  the  covariance 
matrices  of  the  complex  multivariate  normal  populations  with 
Pearson* s  Type  I  distributions.  But  these  approximations  are 
based  upon  empirical  investigations,  whereas  the  investigations 
in  the  present  paper  are  analytic  in  nature.  In  the  real  case, 
Rao  (1951)  gave  a  useful  approximation,  in  terms  of  beta  series, 
for  the  null  distribution  of  certain  power  of  the  multivariate 


beta  matrix. 
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2.  PRELIMINARIES 


In  this  section,  we.-  define  some  notation  and  give  some 
lemmas  which  are  needed  in  the  sequel. 

Let  f(x)  be  a  function  of  the  real  variable  x  defined 
for  x  >  0.  Then,  the  Mellin's  integral  transform  of  f(x) 
is  defined  (e.g.,  see  Titchmarch  (1937))  as 


M{  f  ( •  )  1 1 } 


f ( x ) dx 


(2.1) 


o 

where  t  is  a  complex  variate.  If  f(x)  is  absolutely  contin¬ 
uous  in  (0,1),  then 


f(x) 


i-c+i00  . 

M{ f ( ♦ ) |t}x  dt. 

c-i°° 


(2.2) 


Note  that  when  f(x)  =  (l-x)g~^,  0  <  x  <  1  .f 


M{f ( • ) | t) 


f1xt“1(l  -  x)g_1dx 
0 


r(t)TCfs) 

r(t+g) 


for  Real  (t)  >0  and  Real  (g)  >  0.  Hence 


1 

2tt1 


c+i°° 


C-ioo 


r(t) 

r(t+g) 


x'tdt 


(1  -x)  g“1i 

r(g) 


c  >  0 . 


(2.3) 


The  following  lemma  was  proved  by  Nair  (1940): 


Lemma  2.1.  Let  $(t)  = 


x  p(x)dx  where  p(x)  is  the  density 


of  the  random  variable  x.  If 


<f(t)  =  0(t  v)  (2.4) 

with  Real  (t)  tending  to  infinity,  then  <j)(t)  can  be  expanded 
as 
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<b(t)  =  7  R.  T(  t+a)/T(  t+a+v+i  ) 

i=0  1 


(2.5) 


where  a  is  any  constant. 

Lemma  2.2. 

°o 

Consider  the  series  [  a.x1  which  converges  to  the 

i=l  1 

function  g(x)  in  the  neighborhood  of  x=0  (or  the  asymptotic 


expansion  of  g(x)  when  x=0).  Then 


exp{g(x) }  =  1  +  l  S.x1 

i=l  1 


where 


3 j ' s  satisfy  the  recurrence  relation 


BJ  =  J  “  1  * 


(2.6) 


(2.7) 


We  use  the  following  notations  as  defined  in  James 

(1964).  The  complex  multivariate  gamma  functions  f  (a) 

P 

and  r^(a,K)  are  given  by 

?  (a)  =  ttP^P*1)/2  n  T(a-i+l)  .  (2.8) 

P  i-1 


rp(a,<)  =  rP^P-1)/2 


P 

n  T(a-i+l+k. ) 
i=l  1 


(2.9) 


where  <  =  (k, ,k0,...,k  )  is  a  partition  of  the  integer  k  such 
that  k, >...>k  >0  and  k  =  k, +...+k.  The  transpose  and  con- 

1—  —  p—  ip 

jugate  of  a  complex  matrix  B  are  denoted  by  B'  and  B  respec¬ 
tively.  Also,  let  C<(A)  denote  the  zonal  polynomial  of  a 
Hermitian  matrix  A,  (i.e.,  A  =  A  ) 

Lemma  2.3.  For  any  integer  r,  variate  x  and  Hermitian  positive 
definite  V,  we  have 


lit— 18  lit— 18  II  c— a  0 
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»  xkC  (V) 

l  l 

k=r  k  (k-r) : 

k: 

r 

k=o  k  k: 


-  xr  (tr  V)r  etr  (xV)  (2.10) 


r  r  XKa  (<)C<(V) 

I  2  - - - - -  =  (x2tr  V2- 


x  tr  V)  etr(xV)  (2.11) 


x  a^(<)C  (V)  2  2  r  r+1  r+1 

l  - - - - -  =  {x  2tr  v  ( tr  V)r-xr+1( tr  V)r  1 

r  <  ( k-r) .' 


+  2r  xr+1  tr  V2( tr  V)r  1-rxr(tr  V)r 


+  r( r-1 )xrtr  V2(tr  V)r_2}etr(xV)  (2.12) 


xka2(<)C  (V)  4  9  0  o  o  -3  <? 

I  - - - - -  =  (x  (tr  0^  +  4xJ  tr  V  -2x  tr  V  tr 

o  <  k: 


+  3x2(tr  V)2-4x2  tr  V2 + x  tr  V}  etr(xv) 


(2.13) 


X  3-«}(k)C  (V)  o  o  n  9  9  9 

V  - - - - -  =  {2x  tr  +3xi tr  V)  -3x2tr  V 

0  <  kl 


+  2x  tr  V}  etr(xV) 


( 2 . 14  ) 


where 


in(<)  =  ?  k.(k  -2j) 

1  j=l  J  J 


a2(<)  =  2  V  k.(k2-3jk.+3j2) 
2  j-i  J  i  J 


(2.15) 


The  above  lemma  was  proved  in  Hayakawa  (1972).  We  need  the 
following  in  the  sequel. 
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From  Hayakawa  (1972),  Lemma  1 

(a1(K)+k)C<(V)  =  tr(A3)2CK(V) |v  =  A  (2.16) 

{3a2(K)-2a2(K)+6ka1(<)-6a1(K)+3k2-2k}C<(V) 

=  [8(tr(A3)3)+3(tr(A3)2)2]CK(V)  |y=A  (2.17) 


where  A  is  a  diagonal  matrix  of  eigenvalues  of  V,  3  is  the 
Hermitian  differential  operator  matrix 

3  -  Jg+iSj.  3  =  (3a6),  3r  -  (3*e).  3j  -  ojg) 


and 


3 


a  .1 

- - -  - 


aS  2  3VR  ’  aB  2  3V1 ft 

a8  aB 


where 


V  =  (V  )  =  (VR.+iVIa) ,  6 

aS  aB  aB  a6 


1,  a=B 
0 ,  a^B 


By  definition  of  zonal  polynomial 


l  c  m  =  (tr  V)J 
< 


(2.18) 


and  Eq.  (2.16) 


7  a1(<)C<(V)  =  t r ( A 3 ) 2 ( t r  V)k|y=A  -k(trV)! 


=  k(k-l)  tr  V2( tr  V)k*2-k(tr  V)k. 


(2.19) 
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Multiply  both  sides  of  Eq.  (2.16)  by  a^(c)  and  sum  over  k, 
by  applying  Eq.  (2.19),  we  have 


l  a^(<)6<(V) 
< 


From  Eq.  (2.17) 
l  a2(<)CK(V) 


=  ( 3k2-2k) ( tr  V)k-2k2(k-l)  tr  V2(tr  V)k-2 
+  4  k(k-l ) (k-2)  tr  V3(tr  V)k_3 

+  k( k-1 ) ( k-2 )(k-3)(tr  V2)2(tr  V)k~4  .  (2.20) 

and  Eqs.(2.18)  ,  (2.13),  (2.20),  we  obtain 
=  (3  k2-k) ( tr  A )k  -  3k( k-1 )  tr  V2(tr  V)k~2 


< 

+  2  k(k-l)  (k-2)  tr  V3  (tr  V)k-3  .  (2.21) 


Eqs .  (2.18)  -  (2.21)  are  needed  subsequently.  Here  we  note 

that  Lemma  2.3  follows  by  multiplying  both  sides  of  Eqs. 

k 

(2.18)-  (2.21)  by  and  summing  over  k. 
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3. ASYMPTOTIC  DISTRIBUTION  OF  THE  LIKELIHOOD  RATIO 
STATISTIC  FOR  TESTING  MULTIPLE  INDEPENDENCE 

Let  Zf  =  be  distributed  as  a  complex  multi¬ 

variate  normal  with  mean  vector  'H'  and  covariance  matrix 
Z.  Also,  let  E{ ( Z^-p^ ) (Zj-Lj ) * }  =  ^±y  We  assu™e  that 

Z^  is  of  order  p^  *  1  and  b  =  p^+...+p^#  Consider  the  hypo¬ 
thesis  E  where 


V  Ei-j  =  0  (i«-l . d).  (3.1) 

Let  (Z^ . Z^),  C  j  *  1  > 2  ,  .  .  .  ,N)  .  be  N  independent  obser¬ 

vations  on  Z'  .  Also,  let  A  =  (A,  ),  where 


Aim- 


5i." 


H  m'1 

J  =  1  J 


d). 

(3.2) 


The  likelihood  ratio  statistic  for  testing  Hq  is 


X  =  |  A  |  {  I  |A  J}'1 
i=l  11 


(  3.3) 


1  /  S 

Make  the  transformation  u  =  \  ,  where  s  is  a  constant 

to  be  chosen  to  govern  the  rate  of  convergence.  The  hth 
moment  of  u  is 


E<uh)  -  {  n  (  J  /  ,3.4) 

j  =  l  r(n-j  +  l)  i=l  a=l  F(n+t-a+i) 

S 


By  using  the  Mellin's  inverse  transform,  the  density  of 


u  becomes 
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f(u)  = 


, ,  ,  x  c  +  i00  U  T(n+ j  +  1) 

K(b,d,Pi,a)  f  u.b.!  j.i  s _  dh 


2tt  i 


C-10O 


d 

n 


Pi 

n 


(3.5) 


i=l  a=l 


T(n+^  -  a+  1 ) 


d  b 

where  K(b,d,p  n)  =  IT  ?  F(n-a+l)/  II  r(n-j+l).  Set 

i=l  a  =1  j  =  l 

n=m+5,e  =  c  +  ms  where  n  =  N-l  and  =  —  where  ;  is  also 

s  s 

a  converging  factor  to  be  chosen.  Then  we  have 


e+i°° 


K(b,d,p.,n)  x 

f(u)  =  - - -  usm  1 


u  t  ? ( t )  dt 


(3.6) 


2Tii 


e-ico 


and 


*(t)  =  it  r(|+5-j+i)/  n  n  r(£+6-a  +  i).  (3.7) 

j=l  s  i=l  a=l  s 


By  the  use  of  formula  for  the  asymptotic  expansion  of  gamma 
function,  we  have 


log  T(x+g) 


log  *2~ 


+  (x+g  -  |) 


log  -X  -  X 


CO 


y  (-Dr 

r=l 


Wg> 

r(r+l)xr 


(3.8) 


for  g  bounded  and  Br(g)  is  the  Bernoulli  polynomial  of 
degree  r  .  So , 


OO 

log  <t(  t)  =  log  s;  +  log  t~V  +  l  (3.9) 

r=l  t r 


where 


b(b+l) 

2 


d  p  (p.+l) 

y  — 


V 


(3.10) 
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/  i  n  r  r 
A  -I'1)  s 
r  r(r+l) 


d  pi 


l  l  5  +1( 5-ot+l)-  l  B  ,(5-1+1) 
i=l  =1  i=l  r  1 


(3.11) 


Hence 


r  co  q 

4>( t )  =  sv  t  1 1  +  T 

r=l  t 


(3.12) 


the  coefficient  Q  can  be  obtained  by  the  recursive  equation 


1  r 

Qr  =  r  -  lAl  Qr-Jl’  Q0  = 

C=1 


(3.13) 


Since  $(t)  =  0(t  ' )  ,  we  can  use  Ec .  (2.5)  and  write  ?(t) 
as  follows: 


-v 


Q, 


t”  I  -§}=  [ 

r=l  t  i=0 


r(t+a) 
r ( t+a+v+i ) 


(3.14) 


and  a  is  a  constant  to  be  chosen  to  govern  the  rate  of 
convergence  for  the  resultant  series.  Using  Eq .  (3.8)  to 
expand  the  gamma  function  on  the  right  hand  side  of  Eq.  (3.14), 
we  have 


log 


T(  t  +  a) 

F( t+a+v+i ) 


°o  a  . 

(  V  +  i  )  log  t  +  y  — V- 
j=l  tJ 


(3.15) 


where 


(~DJ 

j(j+D 


CB ,+1( v+a+i) 


Vi(a) 


(3.16) 


Thus 


f(t+a) 

F ( t+a+v+i ) 


t 


-(v+i) 


r  °°  c .  n 

i+  y  -U 
L  jZi  tJ_ 


(3.17) 
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and  can  be  recursively  computed  as 


cij  *  J  J!*AU  ci,j-r  Ci0  -  »• 


(3.18) 


Substituting  Eq.  (3.17)  in  Eq.  (3.14)  and  equating  the  co¬ 
efficient  of  same  powers  of  t ,  FL  is  determined  explicitly 
as 


y 

j*0 


R.  , 


C.  .  .  *  Q. 

i-j,j  i 


1. 


(3.19) 


Now  use  Eqs.(3.14),  (3.12)  in  Eq.  (3.6)  and  note  that  the 
term  by  term  integration  is  valid  since  a  factorial  series 
is  uniformly  convergent  in  a  half-plane  (Doetsch  (1971)). 
Then  we  have 


f(u) 


rr  /  u  r  \  SIT1—  1 

K(b,d,pi,n)  u  s 


I  8. 

j=0  J 


1 

27ti 


e+ic 


-t  r(t+a) 

T ( t+a+v+j ) 


e-i°° 


dt 

(3.20) 


Using  Eq.  (2.3)  for  the  integral,  we  have 


oo 

f(u)  =K(b,d,p  n)  SV  l  R.  usm+a_1  (l-u)v+j_1/r(v+j)  ,0  <  u  <  1  . 

1  j=0  J  “  ~ 

(3.21) 


Thus  the  c.d.f.  of  u  in  terms  of  incomplete  beta  functions 
Ix(-,*)  is 


Prob  (u £ x)  = K(b , d , p . , n)  sv  £  R. 

j=0  J 


Ix(sm+a,v+j  ) 


I*  ( sm+a) 
T(sm+a+v+j  ) 


(3.22) 
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where 


V*’3) 


(a+g)  O- 

r(ct)  r(B)  u 


(l-u)S_1  du. 


(3.23) 


Further  expansion  of  K(b,d,p^,n) 


d  Pi 

n  n  r(n-a+i) 

i=l  a=l 
b 

n  r(n-j+i) 

j-1 


and 


T (sm+a) 
f  (  sm+a+v  +j_) 


gives  us  that 


oo  a; 

log  K(b,d,p  n)  =  log  m J  +  £  — 

1  r=*l  mr 


where 


(  i  \ r  ’  b 
A*  =  (-1) 

r  r(r+l) 


d  Pi 

I  B  (5-i+l)-  l  l  B  .(6-0+1) 
i=l  r  1  i*l  =1  1 


(3.24) 


and 


K(b,d,p.,n)  =  mv  j~l  +  J 

,  r=  1  m 


Q!n 


where 


Qr  =  r  I  i  A£  Q^_  ^ ,  Qq  *  1 


Also 


log 


? ( sm+a) 


A* 


’(sm+a+v+j ) 


rrr  =  log  (ms)  (v+3\+  l  — ir 


(3.25) 


(3.26) 


r=l  m 


where 
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r 

&*-  =  —4—^ -  [Br+1(a+v+j)  -  Br+1(a)]  =  (3.27) 


jr  r 


s  r(r+l) 


and 


r(sm-t-a)  =  ,  v-(v+j) 

r ( sm+a+v+ j ) 


|—  *>  C.  ' 

1+  l  -j£ 

r=l  mr 


(3.28) 


Cjr  -  ?  j/  Aj*t  CJ,r-t  *  ^  Cio  ’  1  •  <3-29> 


Hence,  Eq.  (3.22)  is  of  the  form 

00  Q*  ^  00 

j 


Prob  (u  <  .)  = 


✓  00  <i  00 

1+1-7  I  rh  IY(  sm+a ,  v+j  )  (ms  )-^ 
r=l  nr  J  i=0  J  x 


06  C  . 

l 

r=l  (sm)  j 


(3.30) 


=  I  (sm+a, v  )  +  l  G. 

x  i=l  m1  1 


where 


G. 

i 


=  l  Ri  ,•  I  (sm+a, M+i-j)  l 

• _ A  A  n  _ 


J-0 


z=o 


0  c 


(3.31) 


i-Z 


The  exact  c.d.f.  can  be  calculated  with  the  above  formula 
when  the  sample  size  n  is  small.  A  suitable  choice  of  5 
will  make  m  large  in  order  to  expedite  the  convergence 
of  the  series  in  Eq.  (3.30). 

For  large  sample  size  n,  let  us  now  examine  the  first 
few  terms  of  Eq.  (3.30).  We  know  that 


F 


G1  =  Ix(sm+a,v) 


+  Q*J  +  Rj  Ix(sm+a,v+l)/s  (3.32) 


„  _  x  ,  , ,  :02  ^  Q1  C01  L 

G2  -  Ix(  sm+a ,  v  )  — j  + - +  Q2 

1  s  s  ~ 


+  Rx  Ix(s«+a,v+l) + 


'11 


+  R„  I  (sm+a,v+2)/s' 

Z  X 


(3.33) 


,  C03  Q1  C°2  Q2  C01  * 

G3  =  Ix(sm+a,v)  -3-  +  - 5—  *  -  +  Q3 

L  s  s  s 


G12  ^1G11  G2 
+  Rx  Ix(sm+a,v+l)  J  ^  -  -  +  — 


+  R2  Ix(sm+a, v+2) 


C21  +  Ql 
s3  s2 


+  Ro  I  (sm+a, v+3)/sJ  (3.34) 

o  X 


From  Eqs.(3.24),  (3.25),  we  have 


,  | —  b  d.  Pi  “t 

Q,  =  A*  =  -  k\  l  B«(  6-i+l)  -  l  l  B2(5-a+l)|  (3.35) 

113  Lj  =  i  z  i=l  a=l  — J 


Eqs.  (3.11),  (3.13)  give 


=  A-^  38  —  s A =  sQ^ 


(3.36) 


From  Eq.  (3,19),  we  have 


R1  +  C01  Q1 


(3.37) 


Eq.  (3.16),  (3.18)  give 

C01  *  A01  =  ~2  ^B2(v+a)  “  B2^a^ 


(3.38) 
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2  1 

By  expanding  Bg(x)  =  x  -  x  +  ^,  which  is  the  Bernoulli 
polynomial  and  setting 


b(b+l ) ( 2b+l )  £  p.(pi+l)(2p.+l) 

1  "  6  '  iil  6 


(3.39) 


v1-v 
6o  "  2v 


(3.40) 


aQ  =  ( 1-v  )/2 


(3.41) 


in  Eq.  (3.35)  and  (3.38),  we  obtain 


A*l  =  *1  =  C01  =  0 


(3.42) 


—  0  • 


(3.43) 


Kence  these  results  together  with  Eqs.  (3.11),  (3.13),  (3.16) 
(3.18),  (3.19)  give 


R2  +  C02  ”  *^2 


(3.44) 


C02  A02  S'  [B3(v+aQ)  B3^ao^ 


(3.45) 


Q2  _  A2  =  S“ 


pi 


l  t  B3(S0-“+1>-,I,  B3(5o-i+1> 


i=l  a=l 


i=l 


.  (3.46) 


By  equating  Eq.  (3.45)  and  Eq.  (3.46),  and  expanding 
3  3  2  1 

Bg(x)  =  x  -  ^  x  +  2  x*  we  obtain 


S2  =  (1~v2>v 

2  4  A2 


(3.47) 
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If  we  choose  sq  to  be  the  positive  root  of  Eq.  (3.47),  then 


Now 


Hence 


In  Eq,  (3.19) 


r2  =  0. 


Q2 


0 


R3  =  Q3  ~  C03 


Eqs .  (3.11)  ,  (3.13)  give 


^3  ~  A3 


Sn  = 

o  3 


-  S3  Q* 
o  3 


So 


(3.48) 


(3.49) 


(3.50) 


(3.51) 


(3.52) 


(3.53) 


All  these  identities  induce  that 


Gi  -  g2  =  0 


(3.54) 


R3 

G3  =  T  (Ix(somo+ao'v+3) 


Ix(somo+ao'v)) 


where 


(3.55) 


m  =  n  -  6  . 

o  o 


(3.56) 


i 
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Thus  by  so  choosing  6q,  aQ  and  sQ,  we  have  the  c.d.f.  of 
u  in  asymptotic  form  as 

Prob  (u  <  x)  -  Ix(somo+ao,v) 

+  “I  C Ix(somo+ao’v+3) 
mo  so 

"  Ix(somo+ao’v)]  +  °(mo4)  •  (3.57) 

For  d  =  2,  p1  =  1,  we  have  vQ  =  b-1,  aQ  =  <5q  = 

sq  =  1 ,  and 

Prob  (X  <  x)  =  I  (N-b.b-l)  .  (3.58) 

X 

Table  1  gives  a  comparison  of  the  accuracy  of  the  approx¬ 
imation  by  taking  the  first  term  in  (3.57)  with  the  approxi¬ 
mation  obtained  by  taking  the  first  two  terms  in  (3.57).  In 
the  table,  the  constant  w  is  defined  by  c  =  exp(-w/2) 
where  P[A<c]  =  a  and  the  values  of  w  are  taken  from  the 
tables  of  Krishnaiah,  Lee  and  Chang  (1976).  The  value  of 
P[\<c]  is  denoted  by  a ^  or  according  as  one  term  or  two 

terms  in  (3.57)  are  used.  Also,  Pi=p  for  i=l,2 . d. 
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TABLE  1 

Significance  Level  Associated  with  the  Asymptotic 
Expression  for  the  Likelihood  Ratio  Test  for  Independence 


n 

d=3, 

P=1 ,  a 

=0.05 

d=3  , 

p=2 ,  a=0 

.05 

d=5 , 

p=2 ,  a= 

0.05 

w 

al 

a2 

w 

al 

a2 

r 

w 

al 

a2 

10 

1.459 

.0499 

.0499 

5.143 

.0488 

.0493 

- 

- 

- 

15 

.923 

.0499 

.0499 

2.987 

.0498 

.0499  | 

9.746 

.0454 

.0472 

20 

.675 

.0500 

.  0500 

2.113 

.0499 

6.510 

.0486 

.0493 

30 

.439 

.0502 

.  0502 

_ 1 

|  1.337 

.0495 

.0495 

3.950 

.0496 

.0498 

TABLE  1  (Continued) 

a 

d=3, 

P-1, 

o=0. 10 

d=3, 

P=2 ,  a= 

0.10 

d=5 , 

,  p=2 ,  a 

=0.10 

w 

al 

a2 

w 

al 

a2 

w 

al 

a2 

10 

1.233 

.0998 

.0999 

i  4.679 

i 

.0982 

.0991 

1  — 

- 

- 

15 

.780 

.0999 

.  1000 

|  2.722 

.0995 

.0997  ] 

9.222 

.0928 

.0958 

20 

.570 

.1004 

.1004 

1.926 

.0998 

. 0998  j 

6.167 

.0979 

.0990 

30 

.372 

.0994 

.0994 

1.218 

.0995 

.0995 

3.743 

.0996 

.0999 
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4.  NONNULL  ASYMPTOTIC  DISTRIBUTION  OF  THE  LIKELIHOOD 
RATIO  TEST  STATISTIC  FOR  INDEPENDENCE 
OF  TWO  SETS  OF  VARIABLES 


Consider  the  case  d  =  2,  =  p  and  P2  =  Q  in  Section 

3.  (Z^Zg)  is  distributed  as  complex  multivariate  normal 
with  mean  vector  y’  and  covariance  matrix  Z  where 


11 

Z12 

21 

Z22 

(4.1) 


Also,  A  =  ( -4  jj )  is  as  defined  in  Eq.  (3.2)  and  p  _<  q.  The 
likelihood  ratio  statistic  for  testing  Z ^  =  0  against  the 
alternative  Z f  0  is  given  by 


2  2  2  2  2  ♦ 
where  R  =  diag  and  Yis  are  the 

sample  canonical  correlations  between  and  Z2 . 

Let  n  =  N-l  =  m  +  5  ,  where  5  is  defined  in  Eq .  (3.40), 
00  o 

and  assume 


P 


2 


(4.3) 


2  2  2 

where  W  is  fixed  as  m0_KO°*  Also,  P  =  diag  (p^,...,p  ) 


and 


pi>- 


■>A 


are  the  roots  of  the  characteristic  equation 


"12“22  21 


2_ 

P  £X1 


=  0. 


(4.4) 
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2  2 

Using  the  density  of  ^i»****Yp  given  in  James  (1964)  the 
moment  of  u  =  X  ' S° ,  sQ  as  defined  in  Eq.  (3.41),  be¬ 
comes 


h  9  n  C(P2)  p  C r ( n+k . -i+1 ) V 

E(uh)  =  |l-P2|n  l  l  — -  n  - i - 

k  <  k 


i=l  T(n-i+l)  r(n-q-i+l) 


r ( n-q  +  *  -  i+1) 

_ fo _ 

r(n  +  |  +kri+l) 


(4.5) 


where  (^(M)  is  the  zonal  polynomial  of  Hermitian  matrix 

M,k=  (k^,kg . kp)  is  a  partition  of  k , ^ .  .  .  k^>0  , 

and  k  =  k^+...+kp.  The  non-null  distribution  is  obtained  by 
using  the  inverse  Mellin  transform  and  by  following  the 
argument  from  Eq.  (3.4)  to  Eq.  (3.29).  Substituting  6Q,  aQ 
and  sq  as  in  Eq.  (3.40),  (3.41)  and  (3.47),  we  have 


Prob  (u  <  x)  =  |  I-P2 1 n  l  l  mc 

k  < 


00  Q  (k)  -1 

1  +  T  — - 


r=l  m 


*  ji0  Rj(k)  Ix(somo+ao’v+j+k)(raoso) 


o  -i 
j 


X 


1  + 


00 


l 

r=l 


(4.6) 


Note  that 


(sorao+ao+v+j+k) 

^ ( Somo+a0 )  r(v+j+k) 

s  m  +a  -1  v+i+k-1 

X  U  °  °  °  (1-u)  du  (4.7) 

is  a.  function  of  k,  and  we  have  the  following  expressions 
analogous  to  Eqs.  (3.11),  (3.13),  (3.16),  (3.18),  (3.19), 
(3.24),  (3.25),  (3.27)  and  (3.29) 

l  R  (k)  c.  .  .(k)  =  Q.(k),  Rn(k)  -  1  (4.8) 

j  =0  J  1  J  J  u 


Ix(somo+ao'v+j+k)  = 


Qr(k) 


1 

r 


y  u  (k>  q  (k) 

i=i  *  r“ 


(4.9) 


A  (k) 
r v  ' 


(-l)r  s£  „ 

r(r+l5  Jj  “Vi'VVW’-BrnlV’-1*1”14-!0) 


Clr<k>  '  ?  J1‘Ai*(k)  Ci.r-J(k,:  Ci0<k>  ■  1  (4'U> 

Air(k>  '  Fu+l7  IBr*l<"VkW>  -  W»o»  (4'12) 

Qr(k)  =  r  J  £  A*(k)  Qr-Jl(k)’  Q^k)  =  1  (4.13) 

A*r(k)  ‘  F(7Tl7  jj  CBr+l<6o-i+1>  +  Br+l<So-t>-i*1) 


-  “WW1*1” 


(4.14) 
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CL<k>  *  7  J/A“it<k>  cl,r-i(k>’  C±0<k)  *  1 


(4.15) 


A(r(k>  •  r  i7a7:  CIWVv+k*i)  -  Vl(ao>]  (4'16> 


so  r(r+l) 


and  v  is  defined  in  Eq.  (3.10);  here 


=  pq 


(4.17) 


The  asymptotic  expansion  of  |I  -P 

m 


gives 


i  -  p2 ; n  *  !  i  -  * 

1  1  m 


o  o 


=  exp 


-<W 


tr  w" 


n 


■  r- 1  r  m 


o 


=  exp  (  -  tr  W)  1 1  -  —  (5  +  j  tr  ) 

‘  o 


- o  •;  (-£  tr  W2  +  i  tr  W3) 

o 


2  (io  tr  w  +  l  tr  w2)  >- 


+  0(mQ3) 


(4.18) 


and 


C  (P2 )mk  =  C 
<  o  < 


^mk  =  C  (W) 

J  O  <v  7 


(4.19) 


since  the  zonal  polynomial  C  (M)  is  homogeneous  of  degree 


k.  Furthermore,  the  use  of  formula  for  the  Bernoulli 
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polynomial 


Bn(x+g)  »  r|o(°)  8n.r(x)g’ 


(4.20) 


shows 


A*(k)  =  A*  +  [(25+1)  k+a^K)] 


(4.21) 


c0llk)  ■  r1  -  5i-  U2(a  +v)-l)  k  *  k2] 
o  o 


(4.22) 


Qj(k)  =  Q*-  (o2+5o+i)k  +  \  (  26q+1  )2  k2 


-  J  (26q+1)  ax(<)  +  ^S^nka^K)  +|  a2(<)  a2(<) 

(4.23) 


(a  +v)  j  9 

C*2(k)  =  C*2+—^ - k  +  [(v+aQ+ir  +  (ao+v)]  k(k-l) 

so  2so 


6(a  +v)+8  , 

+  - 2-^ - k(k-l)  ( k-2 )  +  — k  (k-1 )  (k-2)  (k-3)  (4.24) 

12s':  8s„ 

o  o 


Cjx(k) 


ao+v  1 


—  [(2(ao+v)+l)  k  +  k  ]  (4.25) 


where 


A1  =  C01  °’  Q2  +  C02  ~  0 


(4.26) 


as  in  Eq.  (3.42),  (3.50)  ,  by  suitably  chosen  ,  5q ,  sq  in 
Eq.  (3.40),  (3.41),  (3.1?)  and  a^(<)  ,  a^Cx)  are  defined  in 
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Eq.  (2.15).  Also 

R1(k)  =  Q1(k)  -  CQ1(k)  (4.27) 

R2(k)  =  Q2(k)  -  C02(k)  -  (Q1(k)  -  C01(k))C11(k).  (4.28) 

Now  define 

ek 

60(somo+ao’v+j)  =  exp(‘6)  I  kT  rx(somo+ao'v+j+k)  (4'29) 

K 

with  I  (s  m  +a  , v+j+k)  as  defined  in  £q.(4.7)  and 

0  =  tr  W .  (4.30) 

Substitution  of  Eqs .  (4.21)  -  (4  .28 )  in  Eq.  (4.6),  and 

application  of  Eqs.  (2.18)  -  (2.21)  over  <  gives  the 

following  expression  for  the  non-null  asymptotic  distri- 

-2 

bution  of  u  up  to  order  mQ  : 

Prob  (u  <_  x) 

=  35(somo+ao,v) 

+  5T  ’  6e(somo+ao’v)(-'o  tr  W  -  |  tr  W2) 
o 

+  3e(somo+ao’v+1)(25o  -  “5^->  tr  w 

r  a  "HV 

+  3e(somo+ao’v+2)  L-"6°+i^— ^  tr  W  +  tr  "  2s~  (trW)2J 

-s 

+  3a(sorao+ao,v+3)  (tr  W)2  -  |  tr  W*]j 
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[-g2  „ 

Be(somo+ao-v)  [t?  (tr  W)2  •  T  tr  *2  ‘  !  tr  w3 


+  y  tr  W  tr  W2  +  g  (tr  W2)2 


p  a  +v  o  o  6  (a  +V) 

+  Be(somo+ao’v+1)  — )  tr  W"  ( 26o-^T~ 2 — 


’o  °  so 


)  (tr  W); 


a  +v  «n 

"(6o  '  2ip)  tr  w  tr  WJ 


6  46  (a  +v ) 


i —  n  u  **  u  i  a  ~r  v  j  „  v  2(  a  +  V  } 

s  -  -iiTa - 5s— -)  tr 


o  o 


2(a  +v)  ,  9i  o 

( 36o - £ - tr  W2  +  (  ——n(  (v+a+l)2 


O  s  - 

o  o 


+  (v+a  ))  -  — Q— 
oy  7  s 


36  (a  +v)  26 

-°8-°  -^^62  +  ^) 


9 

.  2  .  o 


o  6^  a  +v 
,j  f  O  ,  o 


(tr  W)  +2i  <tr  W)  -(-r  +  2§— )  tr  W  tr  W- 


-  |  (tr  W2)2  +  tr  ff2  (tr  W)2~ ] 

o  J 


7  •>  (an+v)  +  <a+v)  2d  (a+v)  <S  % 

+  MWao’v+3)  -g - 2 - - _j>  tr  W 


r?(a+v)  19^9 

+  (  — ^ -  +  -  45o  J  tr  W  +  |  tr  W 


2(a  +v)  +  3(a  +v)+l  35  (a  +v) 

n  n  a  ' 
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46 

~  +25^+1 
So  ° 


(tr  W)‘ 


a  +v  _ 
0.2 

- K —  + 


35 


)  ( tr  W)3  + 

2so  3so  ° 


55, 


aQ+v+2 


tr  W  tr  W  -  tr  W2  (tr  W2 ) 


+  J  (tr  W2)2 


B3(sorao+ao’v+4) 


3  j.  ~'~o 
2  5o" 


3(a  +v+l) 


2s 


tr  W2  -  2  tr  W3 


r(ao+v)^ + 4(aQ+v)+2  6Q(ao+v+2)  52+l  % 


2s" 


(tr  W)' 


^ao+V  ,  3  35o 


T  +  2i~  -  3s1  (tr  w>  + 


3(a  +v)  +  7 
o 

2s 


-  26 


tr  W  tr  W2  -  tr  W2  (tr  W)2 

so 


+  |  (tr  W2)2  +  -i*  (tr  W)4" 

8s2 

o 


+  VsoVao’v+5) 


o  o  c  a  +v  6 

3  tr  w  +  [  — 2  +  — 2 - 2i~]  (tr  W)3 

6so  2so  ° 


r  6  a  +v+3 
o  o _ 

~  2s 


tr  W  tr  W2  +  tr  W2  (tr  W)' 


-  |  (tr  W2)2  -  -L  (tr  W)41 
4so 


J 
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+  68(som0+ao,v+6)  [-4^-  tr  W2  (tr  W)2+|  (tr  W2)2+-^  (tr  W) 

L  °  8s0 

+  0(nT3)  (4.31) 

When  the  alternative  hypotheses  are  close  to  the  null 
hypotheses,  the  accuracy  obtained  by  using  the  first  two  terms 
in  the  asymptotic  expressions  is  sufficient  for  practical 
purposes.  When  the  alternative  hypotheses  deviate  from  the 
null  hypothesis  significantly,  we  require  higher  order  terms 
to  obtain  very  good  accuracy. 
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5.  NULL  ASYMPTOTIC  DISTRIBUTION  OF  THE  TEST  STATISTIC 
FOR  MULTIPLE  HOMOGENEITY  OF  COVARIANCE  MATRICES 


Let  Z^,...Z  be  independent  complex  p-variate  normal 

variables  with  mean  vectors  p, , . . . ,p  and  covariance  matri- 

-1  ~q 

ces  ,  respectively.  Also,  Let  Z .  .( j  =  l , . . .  ,N.  ) 

i  q  -  i  J  1 

be  independent  observation  on  Z.  .  Let  H  =  n  H  . 

~i  o  j=1  oj 

where  HQj(j=l . d)  is  given  by 


H  .  :Z  *  , ,  =  . . .  =!„* 

°J 


(5.1) 


Qj  =  •  •  • +Qj  .  q*  =  0,  q^  =  q^  and  q^  =  q  . 


(5.2) 


The  modified  likelihood  ratio  statistic  for  testing  Hq  is 


given  by 


n .  d  qj 


n*  -1 

*  I  J  \ 


=  n  |A  /n  |  1  (  n  I  l  A  /n*|  J) 

i=l  1  1  1=1  i=q*  ,+1  1  J 


(5.3) 


where 


n  =N  -1,  n*  =  l  n 

i=qj-l+1 


(5  .4) 


N. 

1  r1 


‘i-  fr  Lin 


!5 .5) 


where  Z  denotes  the  complex  conjugate. 


Let  w  =  A^/n  where  n=  1  and  define  n  =  m  +  6 

i=l  1 
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where  5  is  a  covergent  factor  to  be  determined.  Thus 

3 


ni = (ro+5 


*  rJ  1/S 

)  Y  and  y.  =  l  y.,  u  =  w  '  ,  where  s  is  also  a 

J  i=c^j-l+1 


convergent  factor.  The  h  moment  of  u  is 


E(u  )  = 


n  n 
a=l 


pilY* 

Tot 


/  n 

II  n  s  x 

a=l  a 


p  d  qa  r(y  (m+6)(l+jf-  )-i+l) 


r ( Y*(m-t-  5 )-i+l ) 
r ( Y*  (m+S  )  ( 1*  jg’  )-i  +  l ) 


(5.6) 


The  modified  likelihood  ratio  test  statistic  X  and  its 
moments  are  given  in  Krishnaiah,  Lee  and  Chang  (1976). 

h 

Using  Mellin  inverse  transform  on  E(u  ) ,  and  following  the 
same  lines  as  in  Section  3,  we  have  the  following  expression 
for  the  c.d.f.  of  u: 

00 

Prob  (u  <  x)  =  I  (sm+a,v)  +  l  L-  G.  (5.7) 

x  i=l  m1  1 


where 


Ix(c,d) 


r(c)T(d) 

r(c+d) 


uc  1(l-u)d-1d  u 


o 


(5.8) 


G.  -  l  Rt  .  I  (sm+a ,  v+i-j )  (  \ 

1  -j=0  1_J  x  <1=0  s1 


Q?  Ci-i  i-i 
*  i-  JjJ — -)  (5.9) 


V 


-  <q-d> 


(5 . 10) 


31 


F  l 

r  1=1 


l  A*  C* 

16  Aj£  Cj,r-£ 


C  . 

=  _!£•  c* 

r  ’  LjO 


1 


(5.19) 


where  Br(x)  is  the  Bernoulli  polynomial. 
Now,  let 


1— v 

i  o 

O  2 


n  = 


d 

I 

a-1 


q«  i 

Y* 


Y* 

a 


r  _  n( 2p  -1 ) 
o  ~  6p( q-d ) 


and  let  sq  be  the  positive  root  of 


2  _  (1-Qv 
~24A* 


(5.20) 


(5.21) 


(5 .22) 


(5  .23) 


Then  the  asymptotic  expression  for  c.d.f.  of  u  up  to 
m^,  where  m^  =  n  -  6q  is 

Prob  (u<x)=I(sm  +  a  ,  v ) 
v  —  x  oo  o 

R 

*  4-§  [ysomo  *  vv  +  3>  °-24) 

m0  so 

-  Ix(so%  +  ao’v)]  +  0  (mo4) 


Table  2  gives  a  comparison  of  the  accuracy  of  the  approx¬ 
imation  by  taking  the  first  term  in  (5.24)  with  the  approxima¬ 
tion  obtained  by  taking  the  first  two  terms  in  (5.24).  In  the 
table,  the  constant  w  is  defined  by  c  =  {exp( -w/2) }^/n 
where  P[w<c]  =  a  and  the  values  of  w  are  taken  from  the 
tables  of  Krishnaiah,  Lee  and  Chang  (1976).  The  value  of 
P[w<c]  is  denoted  by  or  ag  according  as  one  term  or  two 

terms  in  (5.24)  are  used.  Also,  q1=...=qd  =  q/d  and  n.^  =  nQ 
for  i-1,2 . q,  n  =  qnQ. 


TABLE  2 

Significance  Level  Associated  with  the  Asymptotic 
Expression  for  the  Likelihood  Ratio  Test  for  the 
Multiple  Homogeneity  of  the  Covariance  Matrices 


n 

o 

q 

d 

P=3, 

Q 

II 

O 

05 

P* 

o 

II 

.05 

w 

ai 

a2 

w 

al 

a2 

6 

3 

1 

37.24 

.0499 

.0500 

67.84 

.0490 

.0493 

6 

5 

1 

64.10 

.0497 

.0499 

118.15 

.0475 

.0486 

10 

3 

1 

33.21 

.0500 

.0500 

56.45 

.0499 

.0500 

10 

5 

1 

57.85 

.0499 

.0499 

| 

100.19 

.0497 

.0499 

20 

3 

1 

i 

30.85 

.0499 

.0500 

50.67 

.0499 

.0499 

20 

5 

1 

54.13 

.0500 

.0500 

90.92 

.0500 

.0500 

10 

6 

3 

46.97 

.0501 

.0500 

81.73 

.0501 

.0500 

10 

6 

2 

58.66 

.0499 

.0499 

102.21 

.0498 

.0499 

20 

6 

3 

43.22 

.0500 

.0500 

72.32 

.  0500 

.0500 

20 

6 

2 

54.49 

.0500 

.0500 

91.76 

.0500 

.0500 

30 

6 

3 

42.12 

.0500 

.0500 

69.73 

.0500 

.0500 

30 

6 

2 

53.26 

.0500 

.0500 

88.85 

.  0500 

.0500 

The  likelihood  ratio  test  statistic  X  given  by  Eq.  (5.3) 


can  be  expressed  as 


\  =  . . . Xd 


(5.25) 


where 


Xj  = 


q3 

IT  |Ai/n.| 


ni 


i=q*+l 

j-1 


K  V»S 


i=q  .+i 
j-i 


"j 


(5.26) 


and  Xj  is  the  likelihood  ratio  statistic  for  testing  Hq ^ . 
An  alternative  procedure  to  test  HQ  is  given  below.  We 
accept  or  reject  ,  according  as 


where 


X  .  >  c 
J  <  <* 


P[Xj  - 


j=l , 


lHo] 


(5.27) 


d 

n 

3  =  1 


P[X  I 

(J 


H  .] 

'  OJ 


(1-a) 


★ 

The  c.d.f.  of  can  be  obtained  by  replacing  d,  q 

and  n  with  1,  q-^  and  n*  respectively  in  Eq.  (5.24);  the  value  of 
is  obtained  by  taking  the  square  root  of  the  right  side 
of  (5.23)  after  replacing  d  and  q  with  1  and  q,  respec- 


I  j  =  2 


l/(Sln*) 

tively.  We  can  similarly  obtain  the  c.d.f.  of  L  J  J 
and  determine  the  values  of  So,  the  value  of 

in  (3.22)  can  be  determined  approximately. 

The  test  procedure  discussed  in  this  section  is  useful 
in  testing  the  hypothesis  that  the  diagonal  blocks  in  the 
covariance  matrix  of  a  complex  multivariate  normal  are  equal 
when  the  off  diagonal  blocks  are  null  matrices. 
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6.  NONNULL  ASYMPTOTIC  DISTRIBUTION  OF  THE  TEST  STATISTIC 
FOR  EQUALITY  OF  TWO  COVARIANCE  MATRICES 


Here  we  shall  consider  the  case  for  d  =  1  and  q  =  2 
in  the  framework  of  Section  5.  The  likelihood  ratio  criterion 
for  testing  the  hypothesis  Hq :  against  the  alternative 

H is  ^ased  on  the  statistic 


X 


npn 

pn,  pnp 
nl  n2 


(6.1) 


where  A=A^+A2 ,  A^  is  as  defined  in  Eq.(5.5),  and  n =  n1+n2 
Let 

1/A 


n  =  mQ  +  6  ■,  n1/n  =  y1>n2/n  =  y2,  w  =  A 


(6.2) 


where  <5q  is  defined  in  Eq.  (5.22)  and  assume  that 

(i)  (I  -  W  is  fixed  as  m0HM:o 

1  o 

(ii)  0  <  lim  <  1;  i=l,2 

1/s 

The  non-null  hx  moment  of  u=w  is 


E(uh)  = 


Y1  Y2 
nl  n2 


es  fph's  n>) 

S  y  l _ Q  ; 

W  fp(n2) 


-n, 


l  l 

k  < 


c  (i-!)"1)  rp<"'k>Winr  +  1>^> 

<  _  _ _ 2 _ 

*  1>’k) 

O 


(6.3) 


(6.4) 


x 
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where  <  =  (k^ , kg , . . . , kp)  is  a  partition  of  k  and  k^^kg^ .  .  .  >_kp_>0  , 
and  sQ  is  the  square  root  of  the  right  side  of  (5.23)  and 
ft  =  ,  while  r  (b)  is  defined  in  Eq .  (2.8),  7  (b,<)  in 

i  p  p 

Eq.  (2.9), 

The  non-null  distribution  is  obtained  by  using  inverse 
Mellin  transformat ion  and  following  the  same  lines  as  in  Sec¬ 
tion  3,  The  c.d.f.  of  u  is  given  by 


c  (I-jT1)  . 

Prob  (u_<  x)  =  |  a  1  _nl  l  l  ~  —  m0 

k  k 


»  Q_(k) 

1  +  l  -2 - 

l  r=1  mn  - 


x  yk  l  Ha(k)  Ix(somo+ao,v+a)(mos0)' 
a=0 


r-  oo  C*  . 

fl  +  I  -T^ 

-  r=l  ml  - 


o 

(6.5) 


where  C*  3  ^  is  defined  in  Eq .  (5.19).  a  *  ~r -  is  defined 

ir  r  o  z 

o  2 

in  Eq.  (5,20)  with  v  *  as  defined  in  Eq .  (5.10).  Also,  we 
have 


i 

.1  Ri_j(k)  ci_jj  =  Q±W,  R0(k)  *  l  (6.6) 

3  0 


r 

QP(k)  =  j  I  Ut(k)  Qr_£(k),  Q0(k)  =  l  (6.7) 

2.  —  1 


and 


A  (k)  * 
r 


(-DrsJ 


7(¥+T)  1=1  Br+l(5o+ki-i+1)  ~ 


Br+l(Y25o-i+1) 


Br+l(YlWi+1) 


yf 


(6.8) 
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Qr(k)  =  r  I  iA*(k)  Qr-Z(k)’  Q0(k)  =  1  (6'9) 


_  (~l)r  £  [  Br+l(Yl5o~1+1)  ,  Br+l(Y2£o  I  +  I) 


Ar('k)  =  rtrTfy  J, 


r 

1 


(6.10) 


-  Br+l<Vkri  +  1> 


From  Eq.  (6.3),  we  have 


--1  _  T  W 

“  "  1  ‘ 

1  o 


hr1!"1-  |i  -  »  f^vV 

1  o 


=  exp 


-  VVV  f  l  tr  W 


i1  ’i 


1-1  ‘"A’1' 


=  exp  (-tr  W)  [1-— (5  tr  W  +  ■—  tr  W2) 

m  o  &  i  -* 

o  1 


_^_{(  otrW2+  1  trw3) 

m  1  3^-i 

o  1 


-  |(5Qtr  W+  ji-trW2)2}] 


+  0(m~  ) 


(6.11) 


Also 


^d-n-1)  Ykmk  .  c  <  *  ) 

1  O 


(6.12) 
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Now  using  Eq .  (4.20)  for  the  formula  of  the  Bernoulli  poly¬ 
nomials,  we  get 


<-l)s0  P  \ 

Al(k)  =  2  ~  JL  B2(5o+ki-i  +  1)" 


B2(Vo~i+1)  B2(YlWi  +  1) 


’1 


(-l)s 


=  A1  +  —2 


-  (1  -  rr-Hk+a.  (<) ) 
Yl 


(6.13) 


A2(k)  =  A2  +  g^jjl  -i)  (-^-2  +  |  +  |  fi1(K)) 

Y1 

+  (l  ~  Y"  )  (35ok+36oa1(<))J 


(6.14) 


Q^(k)  -  Apk)  =  — j  ][ 


B2(Vcri  +  1)  B^(Yo^-i+D 


2V  2  o 


i=l  *-  '1 


-  B2(5  +jt.  -i  +  i)  ! 
*  o  1  1 


(6. 15) 


Q*00 


-  aJ  +  J  ';(2v l)  k  +  ^i(<) 


=  |(A*(k)  Q* ( k )  +  2  A*(k)) 


I- 1  2  2  (26+l)ka,(<) 

*  «5  - 1 5<2Vl)  k  *  - 2 — 4 — — 


12  1  1  1  a9  (  K  )  “1 


(6.16) 


where  ,  A2 ,  Ap  Q2  are  defined  in  Section  5,  and  a1(<), 
a^( k)  are  defined  in  Eq.  (2.15).  By  suitably  chosen  aQ ,  5Q 
and  sQ  as  in  Eq .  (5.20),  (5.22)  and  (5.23)  ,  we  get 
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A1  =  C01  =  °>  «2  +  C02  =  Q2  +  -2§  =  0 


Furthermore,  from  Eq.  (6.6)  -  (6.8),  we  get 


Qx(k)  =  Ax(k) 


Q2(k)  =  J(A1(k)  Q:(k)  +  2  A2(k)) 
Rx(k)  *  Qx(k)  -  CQ1 

R2(k)  =  Q  (k)  -  Cn9  -  (0,  ( k)  -  Cm)C 


(6.17) 


(6.18) 


Ol'^ll 


If  we  now  substitute  these  identities  in  Eq.  (6.5),  apply 
Eqs.  (2.10)  -  (2.14)  for  the  summation  of  zonal  polynomials, 
and  neglect  the  terms  of  higher  order  of  raQ ,  we  obtain  the 
following  expression  for  asymptotic  non-null  distribution. 


Prob(u<x)  =  Ix(somo+ao,v)  +  jjrV  V s0mo+ao , v )-Ix< s0m0+a0 . v+1 )  ] 


+  4  X  aiIx(8omo+VJ+i”1)  +  0(mo3) 


m  i=l 
o 


(6.19) 


where  d^  =  ^  (1  - tr 


5  d,  +  i  d?  +  i  (1  -  4y)  tr  l3 
O  1  2  1  3  .2 


\  (46o  'i±i)  dl"dl“(1_T^  [tr  w3  +  5  (tr  W)2] 
o  1 


1  ,2  ,  ,2  1^1 


a3  =  (5o  "  SI-  +  1)  dl  +2  dl  +  (3  ‘  y7  +  ZT}  tr  W 

O  1  oT^ 


A 
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